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An  Interpolation  Method  with  Weights  and 
Relaxation  Parameters 


Chiew-Lan  Tai,  Brian  A.  Barsky,  and  Kia-Fock  Loe 


Abstract.  This  paper  presents  a new  interpolation  method  that  is  based 
on  blending  a nonuniform  rational  B-spline  (NURB)  curve  with  a singu- 
larly reparametrized  (SR)  linear  spline.  The  resulting  curve  is  called  the 
a-spline.  It  has  weights  and  relaxation  parameters.  Given  the  data  points 
to  be  interpolated,  a NURB  curve  is  obtained  by  using  these  data  points 
as  its  control  points.  The  SR  linear  spline  is  then  determined  by  imposing 
constraints  on  the  a-spline  to  interpolate  the  data  points.  The  a-spline  is 
parametrically  continuous.  It  involves  only  simple  computation,  and  does 
not  require  solving  linear  systems.  This  approach  is  extended  to  produce 
interpolatory  surfaces,  and  can  be  used  as  a modeling  tool  for  deforming 
polygonal  shapes  into  smooth  spline  surface  models. 


§1.  Introduction 

Cubic  interpolating  splines  are  known  to  exhibit  undesirable  oscillations  due 
to  “extraneous”  inflection  points.  This  undesirable  property  motivated  the 
introduction  of  tension  applied  to  interpolating  splines.  Barsky  analyzed  two 
well-known  approaches  for  applying  tension  to  interpolating  curves  [1].  The 
first  approach  is  the  exponential-based  spline  in  tension  [12];  the  second  is  the 
more  efficient  polynomial  alternative,  the  ^-spline  [11].  The  tension  parame- 
ters are  associated  with  data  points  in  these  splines. 

Tension  has  also  been  introduced  as  a shape  parameter  to  non-interpo- 
lating splines.  The  Beta-spline  [2]  incorporates  bias  and  tension  parameters, 
based  on  the  theory  of  geometric  continuity.  Another  spline  technique  that 
has  bias  and  tension  parameters  (and  also  a “continuity”  parameter)  was 
proposed  by  Kochanek  and  Bartels  [8].  Unlike  those  in  the  Beta-spline,  these 
parameters  control  only  the  first  derivatives,  and  the  resulting  curve  is  C 1 and 
interpolatory.  The  weighted  splines  by  Foley  [6]  provide  tension  control  on 
curve  segments  between  interpolation  points,  rather  than  at  the  interpolation 
points. 
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In  addition  to  spline  methods  that  allow  the  user  to  modify  the  shape  pa- 
rameters, a related  research  area  is  the  automatic  determination  of  shape  pa- 
rameters to  produce  curves  and  surfaces  that  satisfy  certain  shape-preserving 
criteria  [5,7].  Another  more  subjective  criterion  used  is  fairness,  or  “visual 
appeal” . 

NURBS  are  also  often  used  to  solve  interpolation  problems.  However,  us- 
ing the  NURB  as  a surface  interpolant  requires  performing  0(mn)  operations 
to  solve  the  0(m  + n)  linear  systems  for  a dataset  of  m x n points.  Although 
these  systems  are  tridiagonal,  they  still  incur  substantial  computations.  More- 
over, the  interpolation  method  is  global;  thus,  changes  to  any  data  point  will 
require  solving  again  all  the  linear  systems. 

The  idea  of  using  singular  blending  to  solve  the  interpolation  problem,  with 
tension  control,  was  proposed  by  Loe  and  Tai  [10].  In  this  paper,  we  use  this 
idea  to  propose  another  interpolation  method  that  has  weights  and  relaxation 
parameters.  A high  relaxation  value  has  the  effect  of  reducing  tautness.  The 
method  is  based  on  blending  a NURB  curve  with  a sequence  of  line  segments, 
or  a linear  spline,  that  is  reparametrized  such  that  each  of  the  line  segments 
has  zero  derivatives  at  both  ends. 

This  work  is  closely  related  to  some  earlier  work  of  Coons  [4].  He  used 
this  blending  technique  to  modify  any  given  piecewise  curve  by  letting  the  line 
segments  be  those  connecting  the  joints  of  the  given  curve,  and  proved  that  if 
the  given  curve  is  Ck,  then  the  blending  function  only  has  to  be  Ck~ 1 for  the 
resulting  curve  to  retain  the  continuity  of  the  original  curve.  More  precisely, 
Coons’  aim  was  not  to  solve  the  interpolation  problem:  for  the  new  curve  to 
be  interpolatory,  the  original  curve  has  to  be  interpolatory. 

§2.  Singular  Blending 

Given  a smooth  piecewise  curve  C(u ),  a tension  control  can  be  introduced  by 
blending  the  smooth  curve  with  a linear  spline  that  approximates  the  curve 
[4,9].  This  linear  spline  can  simply  be  obtained  by  connecting  the  joints  of  the 
smooth  curve,  but  in  general  it  need  not  interpolate  the  smooth  curve.  Let 
L(  u)  denotes  a linear  spline.  The  blending  then  gives 

Q(u)  = (1  — a)C(u)  + aL(u),  it  e [no, un],  0 < a < 1.  (1) 

It  is  easy  to  see  that  as  a increases,  the  contribution  of  L( u)  to  Q(u)  increases; 
thus,  the  resulting  curve  is  more  taut,  simulating  a higher  tension. 

Assuming  that  the  smooth  curve  C(u)  is  C2  (generalization  to  higher 
order  curves  is  straightforward),  then  for  Q(u)  to  retain  this  continuity,  in 
general  L(u)  must  be  at  least  C2.  We  define  L(u)  as 

L{u)  = (1  - s(t))Vj  + s(t)Vj+i,  ue[uj,uj+ 1),  (2) 

where  t = UU~^UJU. , Vj,  j = 0,  are  the  vertices  of  the  linear  spline,  and 
s(t)  is  a monotonically  increasing  function  yet  to  be  defined.  Since  L(u)  must 
be  C2,  it  must  satisfy  the  following  conditions: 

L'(uj)  = L'(uj+ 1)  = 0 and  L"(v,j)  = L"(uj+ 1)  = 0. 
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These  conditions  can  be  satisfied  by  letting  s(t)  be  the  quintic  Hermite  polyno- 
mial 10t3  — 15t4  + 6t5  [3].  We  call  a linear  spline  that  satisfies  these  conditions 
a singularly  reparametrized  (SR)  linear  spline. 

For  the  technique  to  be  useful,  it  should  have  local  tension  control;  that 
is,  the  a in  (1)  must  be  a function  of  u , 

Q(u)  = (1  - a(u))C(u)  + a(u)L(u),  ue[u0,un].  (3) 

The  a{u ) function  must  satisfy  three  criteria:  (1)  interpolate  the  local  parame- 
ters at,  j = 0,  ...,n;  (2)  employ  a local  interpolation  method  so  that  modifying 
a particular  at  will  affect  only  the  neighboring  curve  segments;  (3)  be  at  least 
C2  so  that  the  blended  curve  is  also  C2  (Coons  [4]  showed  that  if  the  SR  line 
segments  are  restricted  to  those  connecting  the  joints  of  the  given  curve,  then 
a(u)  only  need  to  be  C1).  We  observe  that  s(t)  can  be  used  to  define  a(u), 
with  all  three  criteria  satisfied;  hence,  we  define 

a(u)  = (1  - s(t))a*  + s(t)a*+1,  u e [uj,«j+i),  (4) 

where  t = We  then  have  a'(uj)  = a'(uj+i)=a"(uj)=a"(uj+ 1)  = 0. 

A drawback  of  this  definition  of  a(u)  is  that  drastic  changes  in  adjacent  a* 
values  can  cause  the  curve  to  undulate.  An  alternate  method  is  to  estimate 
the  first  and  second  derivatives  at  the  joints  by  some  approximation  method, 
and  obtain  a(u)  by  Hermite  interpolation. 

This  idea  of  singular  blending  can  be  applied  to  many  combinations  of 
curves  and  SR  linear  splines.  Since  NURBS  are  prevalent  in  industry,  we 
let  the  smooth  curve  C(u)  be  a NURB  curve  and  call  the  resulting  curve 
the  a-spline.  If  the  SR  linear  spline  is  the  control  polygon  of  the  NURBS 
curve,  then  the  resulting  a-spline  is  non-interpolatory.  The  non-interpolatory 
a-spline  includes  the  NURB  as  a special  case;  when  all  tension  parameters 
are  zero,  the  a-spline  reduces  to  a NURB.  Some  other  geometric  properties 
inherited  from  NURB  include  convex  hull,  affine  and  projective  invariance, 
and  local  control. 


§3.  The  Interpolating  a-spline  Curve 

The  a-spline  is  non-interpolatory  when  the  SR  linear  spline  is  the  control 
polygon  of  the  NURB  curve.  However,  there  is  no  reason  to  restrict  the  SR 
linear  spline  to  be  the  control  polygon.  In  this  paper,  given  a NURB  curve, 
we  will  determine  a new  SR  linear  spline  such  that  when  they  are  blended, 
the  resulting  a-spline  interpolates  a given  set  of  data  points. 

Let  the  data  points  be  Pj,j  = 0,  ...,n.  The  cubic  NURB  curve  must 
somehow  approximate  the  data  points.  A simple  way  to  achieve  this  is  by 
letting  the  given  data  points  serve  as  the  control  points.  Next,  we  must 
introduce  two  new  control  points  so  that  the  number  of  NURB  curve  segments 
is  n,  where  each  curve  segment  corresponds  to  an  interval  between  two  data 
points.  Since  we  want  the  NURB  curve  to  interpolate  the  endpoints  P o and 
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Pn,  we  let  the  new  points  be  coincident  with  the  two  endpoints,  and  set  the 
first  four  and  last  four  end  knot  values  to  be  equal;  that  is, 


P-i  = P0  and  P„+ 1 = P„, 


and  u_3=ii_2='u_i='ito,  and  un=un+1  =un+2=un+3  in  [/  = {u_3,  ...,un+3}. 
The  interior  knot  values  uq,  ...,  un  are  computed  using  methods  such  as  the 
chord-length  or  centripetal  method.  All  the  examples  in  this  paper  use  the 
centripetal  parametrization.  The  cubic  NURB  curve  sequence  is  given  by 

n+1 

C(u)  = ^2  Rjti(u)Pj,  u g [u0,un], 

j—- 1 


where  Rj,i{u ) are  the  cubic  rational  B-spline  basis  functions  defined  over  the 
knot  vector  U. 

To  determine  the  unknown  vertices  Vj , j = 0, ...,  n,  of  the  SR  linear 
spline,  we  impose  constraints  on  the  resulting  a-spline  Q{u)  to  interpolate  all 
the  data  points:  Pj  = Q(uj),j  = 0,  Substituting  Q(uj)  from  (3), 


Pj  = (1  - a{uj))C(uj)  + a{uj)L{uj)  = (1  - a*)C(Uj)  + a*Vj. 

amP 

Solving  for  Vj,  then  adding  Pj  and  subtracting  ;Q.  ’ yields  equations  that 
require  only  simple  computations: 


Vi 


P3  - (1  - Q*j)C(uj) 


= Pj  + 1—P-(Pj  - C(Uj)). 


Defining  pj  = 


yields 


Vj  = Pj  + Pj  {Pj  - C{uj)),  j = 0, n. 


(5) 


We  call  pj  the  relaxation  parameters.  The  geometric  interpretation  of  (5) 
is  that  C(uj),  Pj,  and  V j are  collinear,  and  the  distance  between  C(iij)  and 
Pj,  and  between  Pj  and  Vj  are  in  the  ratio  1 : pj.  We  know  that  p3  > 0, 
because  a*  < 1.  From  empirical  study,  we  have  found  0 < Pj  < 4 to  be  a 
useful  range.  The  midpoint  value  pj  = 2 yields  visually  appealing  shapes  for 
most  datasets;  thus,  we  use  that  value  as  the  default  relaxation  value.  The 
effect  of  pj  can  be  interpreted  from  noting  that  a*  = and  observing 

the  role  of  a in  (3).  The  value  pj  = 0 corresponds  to  the  maximum  tension 
a*j  = 1;  when  all  pj  = 0,  only  the  SR  linear  spline  contributes  to  the  a-spline 
and  we  have  a linear  interpolant.  When  the  relaxation  value  pj  increases, 
a*j  decreases,  and  the  contribution  of  the  SR  linear  spline  to  the  a-spline 
decreases;  thus,  the  resulting  curve  is  being  relaxed  locally.  By  using  different, 
relaxation  values,  we  can  easily  obtain  rounder  or  sharper  corners  without 
specifying  multiple  knots  or  multiple  control  points.  Fig.  1 shows  the  effect  of 


Fig.  2.  Effect  of  varying  ps:  1,  4;  the  other  pj  are  2. 


Fig.  3.  Effect  of  varying  W5:  2,  0.5;  the  other  Wj  are  unity.  All  pj  are  2. 

applying  the  relaxation  globally.  It  can  be  observed  that  while  the  SR  linear 
spline  found  is  dependent  on  the  global  relaxation  value  specified,  the  cubic 
NURB  curve  remains  fixed.  Fig.  2 shows  the  effect  of  varying  the  relaxation 
locally  (while  the  other  part  of  the  curve  has  the  default  relaxation  value). 
Note  that  the  effect  of  each  pj  is  very  localized;  only  the  nearest  two  curve 
segments  are  affected. 

The  effect  of  the  weight  is  less  obvious.  Each  Wj  affects  four  neighboring 
segments,  i.e.,  [uj- 2,^+2]-  When  wj  decreases,  C(uj ) moves  further  from  Pj 
(assuming  pj  is  fixed);  thus,  from  (5),  V j found  is  further  from  Pj.  Since 
the  weights  are  relative  in  nature,  decreasing  Wj  also  causes  C(v,j-i)  and 
C{uj+ j)  to  move  closer  to  the  edge  Pj-2Pj-i  and  PJ+1PJ+2,  respectively, 
and  to  be  on  that  edge  when  Wj=0.  Hence,  when  Wj  = 0,  from  (5)  again, 
Pj- 2,  C(uj-i),  Pj~  1,  and  Vj-i  are  all  collinear,  and  so  are  Pj+ 2,  C(v,j+ 1), 
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Pj+i , and  V J+\ . This  effect  on  the  SR  linear  spline  is  depicted  in  Fig.  3; 
only  W5  is  varied  here.  To  summarize,  decreasing  Wj  causes  the  interpolating 
a-spline  to  be  rounder  near  Pj,  but  causes  its  segments  between  Pj _ 2 Pj - 1 
and  Pj+jPj+2  to  be  more  taut  near  the  points  Pj-i  and  Pj+\. 

The  a-spline  Q(u ) is  clearly  C 2 continuous  since  a(u),  L(u),  and  C(u) 
are  all  C2.  The  first  and  second  derivatives  at  the  knots  are  as  follows: 

q'm  = (i  - «;)<?'(%) + ot'tomuj)  - cm), 

Q"{uj)  = (1  - a*)C"{u3)  - 2a'(uj)C'(uj)  + a"(uj)(L(uj)  - C(Uj)). 

If  a(u)  is  defined  as  in  (4),  then  Q'{uj)  = (1  — a*j)C'{uj)  and  Q"(uj)  = 
(1  — a*)C"(uj).  That  is,  the  derivatives  at  the  joints  of  the  blended  curve 
Q(u ) are  in  the  same  directions  as  their  counterparts  of  C(u). 

§4.  The  Interpolating  a-spline  Surface 

Analogous  to  the  blending  of  a NURB  curve  with  an  SR  linear  spline,  we 
can  blend  a NURB  surface  with  a network  of  singularly  reparametrized  (SR) 
bilinear  patches.  An  SR  bilinear  patch  is  defined  as  follows: 

L{u,v)  =(1  - s(r))(l  - s(t))V ij  + (1  - s(r))s(t)Viyj+1 

+ «(r)(l  -s(f))Vi+i,j  + s(r)s(t)Vi+i,j+i,  (6) 

where  u e [ui,ui+i),v  € [vj,vj+i),  r = ^7^7,  < = , and  the  functions 

s(-)  are  the  Hermite  polynomials  given  earlier.  That  is,  it  is  parametrized  such 
that  its  first  and  second  order  partial  derivatives  go  to  zero  at  the  boundaries: 

duL(u  i,v)  = duL(ui+i,v)  = dvL(u,Vj)  = dvL(u,Vj+ 1)  = 0, 
dlL{uiyv)  = dlL(ui+ i,v)  = d%L(u,Vj)  = d2vL{u,vj+l)  = 0, 
dlvL(uitv)  = 8lvL(ui+ uv)  = dlvL(u,Vj)  = d2uvL{u,vj+l)  = 0. 

The  a-spline  surface  is  then  given  by 

Q(u,v)  = (1  - a(u,  v))S(u,v)  + a(u,v)L(u,v),  u € [u0,«m];  vG[v0,Vn], 

where  a(u , v)  is  the  blending  function  that  interpolates  the  local  tension  pa- 
rameters ajj,  defined  by  an  equation  similar  to  (6). 

To  find  an  a-spline  surface  interpolating  a given  network  of  data  points 
Pij , i = 0,  ...,m,  j = 0,  ...,n,  we  must  first  find  a NURB  surface,  then  deter- 
mine a network  of  SR  bilinear  patches  to  be  blended  with  the  NURB  surface. 
The  NURB  surface  is  defined  by  simply  letting  the  data  points  be  the  control 
points.  As  in  the  case  of  curves,  we  repeat  the  boundary  vertices  and  let 
the  first  and  last  four  knot  values  be  equal.  Repeating  the  boundary  vertices 
along  the  j-index,  then  repeating  those  along  the  i-index,  we  obtain 

P—i7j  — P®,j  and  fm+rj  1 Pm,j  j j — 0,  ...,7l, 

Pi,— 1 — Pi, 0 and  Pii7l+\  — Pi,m  i — 1, m + 1. 
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Setting  the  end  knot  values  to  be  equal  in  the  knot  vectors  L7={u_3,  nm+3} 

and  V - {v-3,...,vn+3}, 

W— 3 ” U— 2 = W_i  = WO,  W m = W m+1  = Wm+2  = um+3i 

W-3  = V-2  = W-l  = WO,  Wn  = Vn+i  = Vn+2  = Vn+3. 

The  interior  knot  values  can  be  determined  using  any  good  parametrization 
method.  The  NURB  surface  is  then  defined  by 

m+l  n+1 

S(u,v)  = •Ri,4(w)Rj,4 (w)-P.j,  we  [w0,wm];  v e [w0,w„], 

*=— i y=— i 

where  Rit 4(11)  and  Rjti(v)  are  cubic  rational  B-spline  basis  functions  defined 
over  U and  V. 

To  determine  the  vertices  V+j’s  that  define  the  SR  bilinear  patches,  we 
impose  constraints  on  the  a-spline  surface  to  interpolate  the  data  points: 

Pi,j  = Q(v,i,Vj) 

= (1  - aij)s(ui>vj)  + i = 0, m\  j = 0, 

Solving  for  V+j,  we  obtain 

ViJ  = ^ 7+  ~ S(Ui,Vj)). 

ai,j 

1— a*  • 

Defining  pij  — yields 

Vi,j  = Pij  +Pi,j{Pij  ~ S(ui,Vj)),  i = 0 j = 0,...,n. 

This  enables  the  direct  evaluation  of  V{j,  and  avoids  the  necessity  of  solving 
linear  systems  required  by  interpolation  with  NURB  surface. 

§5.  Smoothing  Polygonal  Shapes 

In  addition  to  fitting  a smooth  surface  over  a dataset,  the  proposed  interpo- 
lation method  can  also  be  viewed  as  an  interactive  modeling  tool  for  deform- 
ing polygonal  shapes  (with  an  underlying  rectangular  topology)  into  smooth 
objects.  The  vertices  of  the  polygonal  shape  are  the  data  points  to  be  in- 
terpolated by  a smooth  surface.  With  this  modeling  tool,  the  user  only  has 
to  specify  the  polygonal  vertices,  which  are  fewer  than  the  number  of  control 
points  of  most  spline  schemes.  Manipulating  polygonal  shapes  is  also  simple 
and  easy  for  novice  designers. 

Figures  4 and  5 show  some  modeling  examples,  all  of  which  are  obtained 
from  the  same  input  polygonal  shape,  shown  in  the  top  left  corner  of  Fig.  4. 
In  Fig.  4,  the  relaxation  parameters  are  varied:  globally  in  the  top  row,  and 
locally  in  the  bottom  row  (the  third  row  of  vertices  have  their  relaxation 
parameters  varied  while  the  other  vertices  have  fixed  pij  =1).  In  Fig.  5,  the 
weights  of  the  third  row  of  vertices  are  varied,  while  the  other  weights  are 
fixed  at  unity,  and  the  global  relaxation  is  set  at  p,j  = 1 for  all  i,j. 
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Fig.  4.  Effect  of  global  (top)  and  local  (bottom)  relaxation:  0,  0.25,  0.67,  1.5,  2.3. 


Fig.  5.  Effect  of  local  weight  (third-row  vertices):  2.5,  1,  0.5,  0.25,  0. 


§6.  Homogeneous  Representation  of  a-spline 

It  is  well  known  that  the  rational  B-spline  can  be  viewed  as  the  projection 
of  a polynomial  B-spline  in  homogeneous  space.  This  property  is  important 
because  it  implies  that  all  the  algorithms  for  the  polynomial  B-spline  can  also 
be  applied  to  the  rational  B-spline.  The  rational  B-spline  is  given  by 

C{u)  = YJR]Au)P],  (7) 

where  Pj  = ( xj,yj,Zj ) are  the  control  points  and  , 

Wj  > 0,  are  the  weights,  and  Nj^( u)  denotes  the  cubic  B-spline  basis  func- 
tions. The  polynomial  B-spline  curve  in  the  homogeneous  space,  Ch(u),  whose 
projection  yields  the  rational  B-spline  C(u)  is  given  by 

ch{u)  = (8) 


where  Pj^WjXj , wj  Vj , Wj  Zj , Wj  ) , since  C{u)  = (§t^y,  §t^y).  and  the 

x,  y , z and  w subscripts  denote  the  respective  components  of  Ch(u). 
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We  can  show  that  this  property  is  also  true  for  the  a-spline.  Substituting 
C(u ) from  (7)  into  the  alpha-spline  equation  in  (3),  and  rewriting  it  as  a 
rational  function,  yields 

_ (1  - Q(u))  E NjA(u)wjPj  + a(u)C£,(u)L{u) 

Q[U)  C£(u) 

The  denominator  can  be  expressed  as  (1  — a(u))C%(u)  4-  a(u  )C£(u);  hence> 
the  polynomial  form  of  the  a-spline  in  homogeneous  space  is 

Qh(u)  — (1  - a(u))Ch(u ) + a(u)Lh(u),  u £ [uo,un], 
where  Ch(u)  is  given  in  (8),  and  Lh(u)  is  given  by 

Lh(u)  = C*(«)((  1 - a(t))V*  + s(*)V^+1),  u € [Uj,uj+1), 
where  V1-  = (a;' 1)  is  obtained  from  Vj  = (Xj,y'j,Zj). 

§7.  Conclusion 

We  have  proposed  a new  interpolation  scheme  based  on  blending  a non- 
interpolatory  NURB  curve  (surface)  with  an  SR  linear  spline  (SR  bilinear 
patches).  The  resulting  interpolating  a-spline  inherits  the  continuity,  and 
the  affine  and  projective  invariant  properties  of  the  NURB.  The  method  pro- 
vides weight  and  relaxation  control,  involves  only  simple  computations,  and 
supports  the  modeling  paradigm  of  deforming  polygonal  shapes  into  smooth 
spline  surfaces. 

Several  issues  have  yet  to  be  investigated  for  the  a-spline.  One  example 
is  knot  insertion  which  is  useful  for  shape  refinement  and  rendering.  Another 
issue  is  the  automatic  determination  of  the  parameters  to  satisfy  certain  shape- 
preserving conditions. 
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